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ABSTRACT: 


The  problem  of  developing  practical  suboptimal  filters 
for  nonlinear  systems  is  treated  using  a  different  approach. 
The  developed  filter  (El-F)  is  found  to  fill  in  the  gap 
between  the  Kalman  and  the  Extended  Kalman  filters.  A 
numerical  experiment  to  test  the  performance  of  the  developed 
filter  is  conducted  and  the  results  are  shown. 


I.  Introduction: 


Estimation  problems,  and  filtering  among  them,  are  basically 
concerned  with  extracting  the  best  information  from  inaccurate  obser¬ 
vation  of  signals. 

From  the  control  theory  point  of  view,  the  problem  of  estimating 
the  state  of  dynamical  systems  plays  an  important  role.  Very  often  the 
optimal  control  law  sought  for  a  dynamical  system  is  some  sort  of  a 
feedback  of  its  state.  Take  for  example  the  control  of  a  chemical 
process,  a  nuclear  reactor,  maneuvering  of  a  space  craft,  guidance  and 
navigation  problems,  and  the  problem  of  control  and  suppression  of 
structural  vibrations.  Also,  sometimes,  it  is  of  interest  to  know  the 
state  of  a  dynamic  system.  Take  for  example  the  tracking  of  moving 
objects  like  satellites  in  orbits,  and  enery  missiles.  These  are  just 
a  few  examples  of  the  application  of  this  knowledge. 

Fundamentally,  the  conditional  probability  density  of  the  state 
conditioned  on  available  observations  holds  the  key  for  all  kinds  of 
state  estimators.  The  case  of  a  linear  dynamical  system,  with  measure¬ 
ments  linear  in  the  state  variables,  in  the  presence  of  additive 
Gaussian  noise,  and  under  the  assumption  of  full  knowledge  of  the 
system's  parameters  and  noise  statistics,  has  been  optimally  solved. 

In  that  particular  case,  the  conditional  probability  density  is  Gaus¬ 
sian.  A  Gaussian  density  is  characterized  by  only  two  quantities, 
namely  its  mean  and  covariance.  Therefore,  the  optimal  linear  filter 
has  a  finite  state,  the  conditional  mean  and  thf  conditional  covariance, 
and  is  widely  known  as  the  Kalman  or  the  Kalman-Bucy  filter.  The 
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Kalman  filter  provides  the  minimum  variance  unbiased  estimates.  Also, 
the  filter  structures  Is  linear.  Its  gain  and  covariance  can  be  pro¬ 
cessed  independently  of  the  estimate  even  before  receiving  the  observa¬ 
tions.  These  features  make  the  Kalman  filter  desirable  and  easy  to 
Implement. 

Unlike  the  linear  case,  the  situation  for  nonlinear  systems  Is 
completely  different.  The  conditional  probability  density  is  no  longer 
Gaussian  even  though  the  acting  noise  Is  itself  Gaussian.  In  this  case 
the  evolution  of  the  conditional  probability  density  is  governed  by  a 
stochastic  Integral -partial  differential  equation,  Kushner's  equation, 
or  equivalently  by  an  Infinite  set  of  stochastic  differential  equations 
for  the  moments  of  the  density  function.  Therefore,  the  truly  optimal 
nonlinear  filter  Is  of  Infinite  dimensionality,  and  consequently  is  of 
a  little  practical  Interest.  Therefore,  the  need  for  practical  suboptimal 
filters  is  apparent. 

Inspired  by  Kalman's  results,  a  great  deal  of  research  effort  has 
been  directed  towards  extending  the  linear  results  and  developing 
practical  schemes  for  nonlinear  filters.  Developments  have  relied  on 
two  main  schemes.  One  Is  concerned  with  approximations  of  the  system 
nonl Inearl  ties.  The  other  Is  concerned  with  approximations  of  the  con¬ 
ditional  probability  density  function.  Several  practical  suboptimal 
schemes  have  been  developed  and  a  huge  amount  of  numerical  simulations 
have  been  reported.  A  brief  account  and  discussion  of  these  suboptimal 
filters  Is  given  In  Emara-Shabaik  (1979). 


Still,  the  task  of  theoretical  assessment  of  such  suboptlmal 
schemes  -  In  the  sense  of  providing  a  measure  of  how  far  a  suboptimal 
filter  Is  from  being  a  truly  optimal  -  has  remained  very  hard  to 
achieve.  It  Inherits  the  very  same  practical  difficulty  of  the  optimal 
filter  -  infinite  dimensionality  -  that  one  is  trying  to  escape.  There¬ 
fore,  the  support  of  any  such  schemes  has  to  rely  heavily  on  computer 
simulation  and  for  that  same  reason  not  a  single  scheme  can  be  claimed 
always  superior.  There  are  cases  when  a  particular  filter  has  performed 
better  than  others,  while  there  are  other  cases  where  it  has  not.  The 
final  judgement  is  left  to  experience  and  the  special  case  at  hand. 
Consequently,  the  development  of  a  new  practical  scheme  will  add  to  the 
list  of  contributions. 

The  main  theme  of  this  paper  and  Its  companion,  under  preparation. 

Is  to  consider  the  nonlinear  filtering  problem  from  a  different  approach. 
The  approach  taken  here  is  to  consider  the  problem  as  the  combination 
of  approximating  the  system's  description  and  solving  the  filtering 
problem  for  the  approximate  model.  As  a  result  some  new  schemes  are 
developed.  The  problem  formulation  and  the  proposed  solution  are  given 
next  followed  by  some  numerical  results. 


II .  Problem  Formulation: 

Consider  the  general  nonlinear  dynamical  system  whose  state  x(t) 
evolves  in  time  according  to  the  following  differential  equation, 

dx(t)  *  [A(t)  x(t)  +  f(x(t),t)]  dt  +  Q^t)  dW(t)  (1) 

*(t0)  *  «0  .  t:t0 

where 

x(t)  e  Rn  is  an  'n*  dimensional  state  vector. 

A(t)  is  an  'nxn'  real  matrix. 

f(x(t),t)  is  an  'n'  dimensional  vector  valued  real  function. 
xQ  e  Rn  is  an  'n'  dimensional  Gaussian  random  vector  (GRV)  with 
E  <x0)  *  x0  *  (2) 

and 

Cov(x0,xo)  a  E  { (xQ  -  x0)(x0  -  x0)'}  ■  PQ  +  (3) 

w(t)  e  Rn  is  an  V  dimensional  Wiener  process,  and 
dW(t)  =  W(t+dt)  -  W(t).  Therefore, 

E  {dW(t)}  *  o  for  all  t  >  tQ  (4) 


Cov(dW(t),dW(t))  A  E  { dW ( t )  dW'(t)>  =  (Idt) 

3 

Where  I  is  the  (nxn)  unit  matrix. 

Q^ft)  is  a  real  matrix,  and 

Q(t)  *  Q*s( t)  Q*5  ( t)  is  a  positive  semidefinite  (nxn) 
matrix. 

*  E{ • }  denotes  the  expected  value  of  (•) 
t  Cov(*,*)  denotes  the  covariance  of  (•)• 
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Also,  consider  the  observations  process  dy(t)  to  be  given  by 

dy(t)  -  [C(t)  x(t)  +  h(x(t),t)]  dt  +  R^t)  dv(t)  (6) 

where 

dy(t)  c  Rm  Is  an  'm'  dimensional  observations  vector. 

C(t)  is  an  .'mxn'  real  matrix. 

h(x(t),t)  is  an  'm*  dimensional  vector  valued  real  function. 
v(t)  e  Rm  Is  an  'm'  dimensional  Wiener  process,  and 
dv(t)  *  V(t  +  dt)  -  V(t).  Therefore, 

E  {  dv(t) }  =  o  for  all  t  tQ  (7) 

and 

Cov(dv(t),  dv(t))  £  E  {  dv(t)  dv'(t)  }=  (Idt)  (8) 

R  (t)  is  a  real  matrix,  and 

R(t)  -  R^t)  R**  (t)  is  a  positive  definite  (nxn)  matrix 

We  assume  that  xQ,  w(t),  and  v(t)  are  all  independent  of  each  other 
for  all  values  of  tatQ.  Also,  the  assumption  that  equation  (1)  satisfies 

the  conditions  for  existence  and  uniqueness  of  solution  given  in 

Arnold  (1974),  and  Jazwlnski  (1970)  is  being  made.  This  means  that  our 

dynamical  system  (1)  adults  only  one  solution  x(t),t>tQ  to  be  its  state 

trajectory  in  the  mean  square  sense.  Furthermore,  it  is  assumed  that 

both  f(x(t),t)  and  h(x(t),t)  are  continuous  in  x(t). 

As  it  is  noticed  from  equations  (1),  and  (6),  the  system  structure 
is  considered  to  be  composed  of  two  parts,  a  linear  part  plus  a  non¬ 
linear  part.  Furthermore,  we  assume  that  the  system  behavior  is  dominated 
by  its  linear  paft,  That  is  to  say, 

•8f(x(t),t)|l  <  H  A(t)x(t)||  (9) 
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and 

||h(x(t),t)||<  y  C(t)x<t)||  (10) 

where 

||zj|  Is  the  norm  of  the  vector  z. 

Equations  (1)  and  (6)  along  with  conditions  (9)  and  (10)  can  be  the 
original  system's  description,  what  is  sometimes  referred  to  as  system's 
with  conebounded  nonlinearities.  Also,  it  can  be  a  representation 
obtained  by  linearization  of  a  nonlinear  system,  where  f(x(t),t) 
and  h(x(t),t)  represent  second  and  higher  order  terms.  In  this  case 
conditions  (9)  and  (10)  are  valid  as  long  as  the  system's  state  x(t) 
remains  within  a  small  neighborhood  of  the  nominal  (linearizing) 
trajectory. 

Accordingly,  conditions  (9),  and  (10)  suggest  that  for  a  good  guess  of 
the  system  state  x*(t)  the  following  approximate  equations  for  the 
dynamics  and  observations  can  be  written  as 

dX](t)  ■  [  A(t)  X'j  (t)  +  f(x*(t).t)  ]  dt  ♦  Q^t)  dw(t)  (11) 

dy(t)  -  [C(t)  x,(t)  +  h(x*(t) ,t)  ]  dt  +  R**(t)  dv(t)  Cl 2) 

By  virtue  of  continuity  of  the  nonlinearities  in  x(t),  we  should  note 
the  following.  As  x*(t)  approaches  Xj(t),  the  approximate  description 
given  in  (11),  and  (12)  approaches  the  true  description  in.(l),  and  (6). 
In  fact,  the  following  equation 

dXj(t)  *[A(t)xj(t)  ♦  f(x,Ct),t)  ]  dt  +  fyt)  dw(t), 

*<»„)  ■  *o>  tetD  (,3) 
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and  equation  (1)  have  the  same  solution  both  In  the  mean  square  sense 
and  with  probability  one. 

Thus  follows,  the  filtering  problem  of  the  system  (1),  (6)  can  be 

considered  as  a  unification  of  model  approximation  and  state  estimation 

of  the  approximate  model.  In  other  words,  first  we  approximate  the 

* 

system  description  by  finding  a  suitable  x  (t).  Then,  solve  the  optimal 
filtering  problem  of  the  approximate  model.  The  optimal  filtering  is 
basically  to  seek  the  minimum  mean  square  error  estimate  of  the  state 
x(t)  based  on  the  available  observations,  Yt=[y(s),  tQss>t]. 

Generally,  according  to  theorem  (6.6)  of  Jazwinski  (1970),  pp.  184; 
and  its  specialization  to  linear  systems;  theorem  (7.3)  pp.  219  of  the 
same  reference,  the  optimal  filter  Imitates  the  dynamics  of  the  system 
and  Is  linearly  driven  by  the  net  observations.  Therefore,  guided  by 
these  results,  we  will  seek  the  optimal  filter  for  the  system  in  (11)  and 
(12)  as  a  linear  dynamic  system  driven  linearly  by  the  net  observations. 
The  optimality  of  the  filter  is  In  the  sense  of  achieving  minimum  mean 
square  error. 

so,  if  we  define  the  estimation  error  e^(t)  as 

ej(t)Axj.(t)  -  Xj(t)  (14) 

and  the  covariance  matrix  P(t)  as 

P(t)$  E{(e.j (t)  -  e(t))(ej (t)  -  ^(t))'}  (15) 

Where  x,(t)  Is  an  estimate  of  x^(t)  based  on  Y^,  and 

e-j(t)  -  E  (e.|(t)}  (16) 

then, 

Jfe^t))  -  tr  (E{e1(t)e'1(t)}) 

-tr(P(t))  +  trfeJt)  ef(t)) 

Is  to  be  minimized. 
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Proposed  Solution,  Derivation  of  the  (E-l)  Filter: 


According  to  our  approximation  the  dynamic  measurements  can  be 
written  as  follows. 


dX| (t)  =  [A(t)x1 (t)  +  f(x*(t),  t)]  dt  +  QH(t)  dW(t) 

dy-j (t)  =  dy (t)  -  h(x*(t),  t)  dt  =  C(t)  X^t)  dt 
+  Ris(t)  dV(t) 

Then  a  linear  filter  structure  can  be  sought  as 

d*j  (t)  =  [B(t)  X'j  (t)  +  f (x*(t)  ,t )]  dt  +  K(t)  dy^t) 


where 


B(t)  is  an  'nxn'  matrix 
K(t)  is  an  'nxm'  gain  matrix 


B(t),  and  K(t)  are  to  be  chosen  to  provide  a  minimum  variance  unbiased 
estimate.  By  definition,  the  estimation  error  is 


e^t)  =  Xj(t)  -  (t) 


(23 


Therefore , 


dej(t)  =  dx^j  (t)  -  dxj(t)  (24) 

And  by  substitution  of  (20).  (21).  and  (22)  in  (.24)  above  we  get, 

de^(t)  =  A(t)x-j(t)dt  +  Q*s(t)  dw(t)  -  B(t)  x-j  ( t)  dt 

-  K(t)  C(t)  x-|(t)dt  -  K(t)  R^t)  dv(t) 

=  [A(t)  -  B(t)  -  K(t)  C(t)]  x] (t)dt 

+  B(t)  e-j  ( t) dt  +  QJs(t)  dw(t) 

-  K(t)  R**(t)  dv(t)  (25) 

It  is  desirable  to  have  the  estimation  error  independent  of  the  state. 

In  this  case  large  state  variables  can  be  estimated  as  accurate  as  small 
state  variables.  Therefore,  we  may  choose 

B(t)  =  A(t)  -  K(t)  C(t)  (26) 

Therefore,  the  dependence  of  the  estimation  error  on  the  state  is 
el  iminated. 

Also,  by  choosing 

J(t0)  *  \  (27) 

we  have 

de,  (t)  »  |A(t)  -  K(t)  C(t)J  e^ (t)dt  +  Q^t)  dw(t) 

-  K(t)  R^t)  dv(t) ,  e,(t  )  «  x  -  x_ 


(26) 


de^t)  -  [A<t)  -  K(t)  C(t)]  e^(t)dt,  e^)  =  o 


(29) 


from  which  follows  that 

ej(t)  =  o  for  all  t  l  tQ  (30) 

Hence,  the  estimate  is  unbiased.  Next,  we  seek 

the  filter  <jain  matrix  K(t)  that  provides  the  minimum  variance  estimate. 
By  definition,  the  covariance  matrix  P(t)  is 


P(t)  A  E  {(ei(t)  -  e^tMe^t)  -  e^(t))'}  (31) 

and  due  to  (30) 

P(t)  =  E  {e^t)  ej(t)}  .  (32) 

Straight  forward  manipulations  show  that  P(t)  is  given  by  the  following 
differential  equation, 


dP(t)  »  ([Aft)  -  K(t)  C(t)]P(t)  +  P(t)[A(t)  -  K(t)  C(t)]' 
+  Q(t)  +  K(t)  R(t)K'(t))  dt 


(33) 


As  the  matrix  differential  equation  for  the  covariance  matrix  P(t). 
The  initial  condition  for  (33)  is  given  by 


HtJ  -  E  «X0  -  *0)(x0  -  x0)>  *  P, 


Now,  the  optimization  problem  involving  the  choice  of  the  gain  matrix 
K(t)  can  be  stated  as  follows. 


tr(P(t) ) 


t0<s<t 


Subject  to  conditions  as  given  by  equations  (33)  and  (34)  which  is 


the  same  as 


min  tr(Pft  +  /dP(t)  )  5  tr(P  )  +  min  tr(  /dP(t)  ) 
K(s)  l°  {  °  K(s)  \ 

_ A.  1.  *  C  C  »» 


vsst 


v$st 


Therefore,  we  seek  K(s),  t  ss<t  that  minimizes 


tr( /dP(t)  )  •  ytr(dP(t)  ) 
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substituting  for  dP(t)  in  equation  (35)  we  get 


./£(dP(t))  *y*t|tr([A(t)  -  K(t)  C(t)]P(t)) 


o 

+ 


tr(p(t)[A{t)  -  K(t)  C(t)]') 
tr(K(t)  R(t)  K(t))  +  tr(Q(t)))dt 


(36) 


The  integrand  in  (36)  is  a  convex  quadratic  in  K(t).  According  to  the 
theory  of  calculus  of  variations,  the  minimizing  K(s),  tQ«s<t  is  given 
as  the  solution  of  the  Euler's  equation  which  reduces  to  a  simple 
algebraic  equation  In  the  present  case,  namely 

alcftT  tr(dP(t))  *  o  (37) 


Using  the  concept  of  gradient  matrices  and  the  formulae  developed  in 
Athans,  et.al .(1965),  we  get 


sicftT  tr(K(t>  C(t)  ’’<t»  *  »*(«>  C'Ct) 

(38) 

sicfty  tr(P(t)  C'(t)  K'(t))  ■  P(t)  C'(t) 

(39) 

and 


jj^ytr(K(t)  R(t)  K'(t))  -  2K(t)  R(t) 


(40) 
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Therefore  the  optimal  gain  is  given  by 


-2P(t)C'(t)  ♦  2K(t)  R(t)  *  o 

l.e. , 

K(t)  =  P(t)C'(t)  R”1 (t) 


(41) 


(42) 


Substituting  this  result  into  equation  (33),  the  covariance  matrix 
satisfies  the  following  matrix  Riccati  differential  equation. 


dP(t)  =  [A(t)P(t)  +  P(t)  A'(t)  -  P(t)C'(t)R~1(t)C(t)P(t) 
♦  Q(t)]  dt,  P(tQ)  =  PQ 


To  summarize,  the  solution  for  the  filtering  problem  of  the  approximate 
model  is  given  by 


dXj(t)  *  [A(t)  x,(t)  +f(x*(t),t)]  dt  +  K(t)  [dy(t) 

-  C(t)  (t)dt  -  h(x*(t),t)dt] 

W  *  *0 
K(t)  -  P(t)C'(t)R"1 (t) 

dP(t)  -  [A(t)P(t)  ♦  P(t)  A'(t)  -  PUK'UJR’VWOfU) 
•H)(t)]  dt 

«  P0 
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According  to  the  argument  following  equations  (11)  and  (12),  x  (t) 

Is  required  to  provide  the  optimal  solution  of  the  following  minimization 
problem. 

nix  J(x*(t)>  -  Ev  {(x,(t)  -  x*(t))'(x,(t)  -  x*(t))}  <«) 

x*(t)  * 

then  for  every  t  >  tQ  setting  3J(x*(t))/3x  (t)  *  o  we  get 

x*(t)  «  Ey^  {  x, (t) J  -  x, (t)  (46) 

Therefore,  combining  the  results  of  equations  (44),  and  (46)  we  get 
the  first  of  the  developed  filters,  to  be  denoted  as  the  (El-F)  filter 
namely, 

d^(t)  *  [ A( ‘t )x-|  (t)  +f(J1(t),t)]dt  +  K(t)  [dy(t) 

-C(t)  5,(0  dt  -hCx, (t) ,t)  dt]  .  x,^)  -  x0  (47) 

K(t)  *  P(t)C'(t)R’1(t)  (48) 

dP(t)  *  [jA( t )P( t)  +  P( t)A"(t) .  -  P(t)C'(t)R*1(t)C(t)P(t)  +Q(t)]  dt 

P(t0)  -  P0  (49) 

It  Is  straightforward  to  recognize  that  In  case  of  a  linear  system, 
l.e.  f(x(t),t)  and  h(x(t),t)  are  Identically  zero  or  only  functions  of 
time,  equations  (47),  (48)  and  (49)  reduce  to  the  well  known  Kalman 
filter. 
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Although  the  (El-F)  Is  developed  using  a  different  approach,  It  bears 
a  close  relationship  with  the  extended  Kalman  filter  (EKF)  given  in  Jazwinski  (197 
The  equations  for  the  state  estimate  of  both  the  (El-F)  and  the  (EKF)  have 
the  same  structure.  While  the  equations  for  the  gain  and  covariance  of 
the  (El-F)  are  different  from  those  for  the  (EFK),  they  are  identical 
to  those  of  the  Kalman  filter  (KF).  Therefore,  unlike  the  (EKF),  the 
gain  and  covariance  for  the  (EL-F)  can  be  processed  off  line  and  prior 
to  receiving  the  observations.  This  is  due  to  the  fact  that,  the  matrices 
A(t)  and  C(t)  in  the  (El-F)  are  different  from  the  corresponding  matrices 
A(x(t),t)  and  cf(x(t),t)  In  the  (EKF),  they  are  no  longer  estimate  dependent. 

The  two  sets  of  matrices  are  related  as  follows. 


ACx(t).t)  »  A(t)  ♦  I  X(t)  =  x(t) 


?<i(t).t)  -  C(t)  ♦  I  x(t)  •  x(t)  (51) 

Therefore,  the  (El-F)  has  the  gain  and  covariance  computational 
facility  enjoyed  by  the  linear  filter,  and  moreover.  It  Is  of  higher 
sophistication  since  it  accounts  for  otherwise  neglected  nonlinearities. 


Therefore,  the  El-F  will  be  of  an  advantage  over  the  EKF  when  on 
line  computations  of  the  gain  and  covariance  are  not  affordable  due  to 
capacity  limitations  of  on  line  computers.  This  is  usually  the  case  of 
airborn  and  spaceborn  computers. 


Furthermore,  while  the  (EKF)  has  to  be  strictly  Interpreted  In  the 

A 

Ito  sense,  Emera -Shaba Ik  (1980),  It  Is  not  the  case  with  the  (El-F). 

This  Is  so  because  the  gain  K(t)  as  given  by  equation  (48)  Is  not  estimate 
dependent.  On  the  other  hand.  If  equations  (1)  and  (6)  are  obtained 
through  linearization  of  some  nonlinear  system,  where  the  system  Is  being 
continuously  rellnearlzed  around  the  most  recent  available  estimate  then 
the  (El-F)  and  (EKF)  are  Identically  the  same.  So,  in  a  sense  the  (El-F) 
provides  the  missing  link  between  the  Kalman  and  extended  Kalman  filters, 
and  this  credits  the  new  approach  of  looking  at  nonlinear  filtering. 


IV.  Numerical  Experiment: 


The  Van  der  Pol  oscillator: 

The  Van  der  Pol  oscillator  Is  characterized  by  the  following 
differential  equation,  Cunningham  (1958). 


x(t)  -  ex(t)(1  -  x2(t))  +  x(t)  *  o 


(52) 


which  describes  a  dynamical  system  with  state  dependent  damping  coefficient 

p 

equals  -  e(l-x  (t))  where  e  Is  a  positive  parameter.  The  damping  in 
the  system  goes  from  negative  to  zero  to  positive  values  as  the  value 
of  x2(t)  changes  from  less  than  to  greater  than  unity.  The  oscillator's 
response  is  characterized  by  a  limit  cycle  In  the  x(t),  x(t)  plane  (the 
phase  plane).  The  limit  cycle  approaches  a  circular  shape  as  c  becomes 
very  small.  It  has  a  maximum  value  for  x(t)  equals  2.0  irrespective  of 
the  value  of  e.  This  type  of  oscillations  occur  in  electronic  tubes 
which  exhibit  also  what  is  known  as  thermal  noise.  Denoting  x(t)  as 
*1 (t) ,  and  x(t)  as  x2(t),  equation  (52)  can  be  rewritten  In  a  state 
space  formulation.  Also,  considering  the  existence  of  some  noise 
forcing  on  the  system,  we  get. the  following  representation  for  the 
Van  der  Pol  oscillator. 


i 
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Also  suppose  that  the  following  measurement  Is  taken 

4y(t)  -  [xl(t)  +  x*(t>]  dt  +  Rh  dv(t)  (5< 

In  (53)  and  (54)  above  [Wj(t)  Wg(t)]T  Is  considered  to  be  a  two 
dimensional  Wiener  process.  Also.  V(t)  Is  a  one  dimensional  Wiener 
process.  R  Is  a  positive  nonzero  real  value,  and  Q  Is  a  (2x2)  matrix 
The  following  values  for  noise  statistics  are  considered. 


Also  e  Is  taken  to  be  0.2 

In  the  figures,  the  following  symbols  are  used. 

XI  =  the  1—  state ,  1-1,2 

XIK  s  the  estimate  of  the  1—  state  provided  by  the  (K-F) 

XIE  s  the  estimate  of  the  1—  state  provided  by  the  (El-F) 

XIEK  =  the  estimate  of  the  1—  state  provided  by  the  (EKF) 

In  both  cases;  as  Indicated  by  figures  1,  2,  3,  and  4,  both  the  (El-F) 

and  (EKF)  provide  very  accurate  tracking  of  the  system's  states  while 
the  (KF)  provides  crude  estimates. 


V.  Conclusions: 

A  new  approach  for  nonlinear  filtering  Is  developed.  Basically, 

It  consists  of  a  model  approximation  technique  combined  with  optimal 
filtering  of  the  approximate  model.  The  resulting  nonlinear  filter 
(El -F)  has  a  structure  that  fits  into  the  gap  between  the  Kalman  and 
the  Extended  Kalman  filters.  On  one  hand  It  enjoys  the  same  computational 
facility  for  the  gain  and  covariance  enjoyed  by  the  Kalman  filter  (KF). 
While  on  the  other  hand  It  provides  estimates  on  the  same  level  of 
accuracy  as  provided  by  the  ectended  Kalman  filter  (EKF). 
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Figure  2.  Second  state  and  estimates  by  Kalman,  El,  Extended 
Kalman  filters. 
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Figure  3.  First  state  and  estimates  by  Kalman,  El,  Extended 
Kalman  filters. 
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Figure  4.  Second  state  and  estimates  by  Kalman,  El,  Extended 
Kalman  filters. 
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